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bJQ. Abstract 

■ The gravitational collapse of a cylindrically distributed perfect fluid is studied. 

^ ■ We assume that the collapsing speed of the fluid is very large and investigate such 

a situation using recently proposed high-speed approximation scheme. We show that 
if the value of the pressure divided by the energy density is bounded below by some 
positive value, the high-speed collapse is necessarily decelerated by pressure repulsion 
so significantly that this approximation scheme becomes inapplicable. However, even 
in the case of a mono-atomic ideal gas, which is a typical example of the bounded case, 
an arbitrarily large tidal force for freely falling observers can be realized before the 
high-speed approximation breaks down if the initial collapsing velocity is set to a very 
large value. By contrast, in case that the equation of state is sufficiently soft, we find 
that the high-speed collapse is maintained until a spacetime singularity forms. 
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§1. Introduction 



Gravitational collapse is an important topics in general relativity. With recent obser- 
vational studies of gravitational waves using laser interferometric detectors such as LIGO, 
TAMA, VIRGO and GEO,^)'^) the study of gravitational collapse as a source of gravita- 
tional radiation was gained its importance. Almost all theoretical studies on generation 
mechanisms of gravitational emissions have been proceeded under the assumption that the 
spacetime singularities are hidden inside black holes. However, if the spacetime singular- 
ities are not enclosed by event horizons, the situation could be completely different from 
the case with horizons. Such a spacetime singularity is called a globally naked singularity. 
Nakamura, Shibata and one of the present authors have conjectured that a large spacetime 
curvature in the neighborhood of a globally naked singularity can propagate away to infinity 
in the form of gravitational radiation due to the lack of an event horizon and therefore the 
mass of the naked singularity is lost through large gravitational emissions. 

Here we should note the cosmic censorship conjecture.^^ The statement of the cosmic 
censorship is, roughly speaking, that a spacetime singularity formed as the result of physically 
reasonable initial data cannot be naked. In connection with this matter, recently, Harada 
and one of the present authors proposed the alternative concept of the spacetime singularity, 
called a 'spacetime border'. ^'^^ There is a cut-off energy scale A above which general relativity 
cannot describe any physical process; for example, A could be the Planck scale Epi (~ 
lO^^GeV), or it might be a scale on the order of several TeV in the brane world scenario. ^^^"^^^ 
Thus, spacetime regions of energy scales higher than A should be regarded as singularities 
for general relativity. This viewpoint leads to the definition of a spacetime border as follows. 
A spacetime border is defined as a region of the spacetime which satisfies the inequality 

infF>^, (M) 

where the curvature strength F is given for instance by 

F :— max(|i?"a|, \R°'''Rab\^^'^ , \R"'^'^'''Rabcd\^^'^, \RiJ.up<j\)- (1"2) 

Here we have extended the original definition of F to a slightly more general one. In this 

case, F includes components Rn^pa of the Riemann tensor with respect to the tetrad basis 

transported parallelly along timelike curves of bounded acceleration, i.e., the tidal force 

experienced by physically reasonable observers. In a practical sense, the visible border is a 

useful concept as an alternative to the naked singularity. 

*) The visible spacetime singularity is called the naked singularity. The globally naked singularity is a 
special case of the naked singularity. It is defined as the naked singularity visible from infinity. 
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In the investigation of naked singularity formation, mainly spherically symmetric systems 
have been studied/^) because they are simple but possess rich physical content. However, in 
a sense, such a system is too simple, as there is no degree of gravitational radiation. There- 
fore, in order to study the generation mechanism of gravitational radiation, we have to add 
non-spherical perturbations to this system^-*"^^ or consider non- spherically symmetric space- 
times. In this sense, cylindrically symmetric gravitational collapse has significant physical 
meaning, because there is a degree of gravitational radiation, and, further, the spacetime 
singularity in this system is naked. ^^^'^^^ There are a few numerical studies of gravitational 
wave generation through cylindrical gravitational collapse. ^^^^^^^ Although there has been 
debate about their results, the present authors have revealed with an analytic study that 
there is no inconsistency.^^-* In any case, these investigations give strong evidence for the 
generation of a large amount of gravitational radiation through cylindrical gravitational col- 
lapse. There are also several studies concerning cylindrical gravitational collapse, and these 
have improved our understanding of non-spherical relativistic dynamics. ^^^"^^^ 

In this paper, we investigate the gravitational collapse of a cylindrical, thick shell com- 
posed of a perfect fluid by generalizing a high-speed approximation for a cylindrical thick 
dust shelP^^ to the case of non-vanishing pressure. In this approximation scheme, the col- 
lapsing speed is assumed to be almost equal to the speed of light. Therefore we treat the 
deviation of the 4-velocity of the perfect fluid from null as a small perturbation. By using 
this approximation scheme, we study the effect of the pressure in the case of a very large 
collapsing velocity. 

This paper is organized as follows. In §2, we briefly review the formulation of a space- 
time with whole-cylinder symmetry. Then in §3, we present the high-speed approximation 
scheme for a general perfect fluid. In §4, we discuss the effect of pressure and determine 
whether the pressure prevents high-speed collapse. Because an ideal mono-atomic gas be- 
haves like radiation in the high energy limit, in §5 we solve the Euler equation with a linear 
equation of state including the radiation case, and show that if the initial collapsing velocity 
is chosen appropriately, then an arbitrarily strong tidal force for freely falling observers is 
realized in the neighborhood of the symmetric axis of the cylinder before the high-speed 
approximation becomes inapplicable. Finally, §6 is devoted to a summary and discussion. 

In this paper, we adopt c — 1 unit and basically follow the convention and notation in 
the textbook by Wald.=^°) 
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§2. Cylindrically symmetric perfect fluid system 



We focus on a spacetime of whole-cylinder symmetry. Such a system possesses the line 
element^^^ 

ds^ = e2(T-V') (-^^2 ^ ^^2^ ^ g2V-^^2 ^ e-2V^2^^2^ (2-1) 

where the metric variables 7, ijj and R are functions of t and r. Then, the Einstein equations 
are 

- SttGv^ (i?'r*t + i?r^) |, (2-2) 

- 87rGx/=^ (rT\ + i^'T'-i) 1 , (2-3) 



7 - 7" = z^'^ - - (2-4) 
R-R" = -SnGy/^ {T\ + T^) , (2-5) 

+ f - r - f V'' = + - T\ + r%) , (2-6) 

where the dot represents differentiation with respect to t and the prime that with respect to 
r. 

In this paper, we consider a perfect fluid. The stress-energy tensor is written 

T^,v = (p + p)Uf,u^ + pg^,,, (2-7) 

where g^^, is the metric tensor, p is the energy density, p is the pressure, and ■u'^ is the 4- 
velocity of the fluid element. Due to the assumption of the whole- cylinder symmetry, the 
components of the 4-velocity u'^ are written 



u 



^^u'(l,-l + V, 0, 0), (2-8) 



where V should be positive so that u'^ is timelike. With the normalization u'^Ujj, = —1, u is 
expressed as 



e 



u* = -^^^^=. (2-9) 
^JV (2 - V) 



Here we introduce the new variables D and P defined by 



" (2 - V) ~ V{2-V) ' ^'-'"^ 
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where g is the determinant of the metric tensor g^y. The components of the stress-energy 
tensor T'^p are then expressed as 

^T\ = e^-'^ {-D + V {2-V)P} , (2-12) 

= e^-^D{l -V) = -^n, (2-13) 

yf^T\ = e^-'^ {{1 - Vf D + V {2 - V) P] , (2-14) 

= v^T% = e^-^y(2 - y)P, (2-15) 

and the other components vanish. 

The equation of motion V^T"^ = leads to 

duD = -\{DV)' + \ {V{2 - V)Py-r ^(1 - V) [2d^{^P - 7) - V{i^ - 7)} 

+ lpV{2-V){^-^-R/R), (2-16) 

Da^V^ = (1 - V^)^^^ + ^ - V^)L> - V{2 - V)P}' 

- ^ [2duii, -^)-V{^-^)] + ^PV{2 - V){y -^' + R'/R), (2-17) 

where u = t — r is the retarded time and du is the partial derivative of u with fixed advanced 
time V — t + r. The first equation comes from the ^-component, while the second one comes 
from the r-component. The z- and (^-components are trivial. 

To estimate the energy flux of gravitational radiation, we adopt the C-energy E and 
its flux vector proposed by Thornc.^^^ The C-energy E = E{t,r) is the energy per unit 
coordinate length along the ^-direction within a radius r at time t. It is defined by 

£; = l|l + e-2T(i?2_^/2^| (2-18) 

The energy flux vector associated with C-energy is defined by 

v/=iJ''=f¥.-f.0.oV (2.19) 



\ (9r ' dt 

By definition, is divergence free. Using the equations of motion for the metric variables, 
we obtain the expression for the C-energy flux vector as 

v^J* = ^^^RR'{^P'' + V'") - 2i?W' - 8nGy^{R'T\ + RT\) | , (2-20) 

= + V'") - 2i?i?W - SnGy/^iR'T't - RT\)^, (2-21) 

while the other components vanish. 



5 



§3. High-speed approximation for a perfect fluid 



Let us consider the ingoing null limit of the cylindrical perfect fluid. Using D and P, the 
stress-energy tensor is written 

7;.. = — ^ {Dk^K + ^^^-"^^¥(2 - V)Pg^J\ , (3-1) 

where 

k^ = {l,-l + V, 0, 0). (3-2) 

The timelike vector k^^ becomes an ingoing null vector in the limit y — > 0+. Hence, in this 
limit with D and P flxed, the stress-energy tensor is identical to that of collapsing null dust, 



Tfii/ * — k^ky, (3'3) 



where 

0, 0). (34) 

This implies that in the case of a very large collapsing velocity, i.e., for < F ^ 1, the 
perfect fluid system is approximated well by a null dust system. For this reason, we treat the 
deviation V of the 4-velocity from null as a perturbation and perform a linear perturbation 
analyses. 

In the case of collapsing null dust, the solution is easily obtained as 

= 0, (3-5) 
7 = 1b{v), (3-6) 
R = r, (3-7) 

SuGDe'' = (3-8) 
dv 

where 7_B(f ) is an arbitrary function of the advanced time v. This solution was obtained by 
Morgan^^^ and was studied subsequently by Letelier and Wang^^^ and Nolan^^^ in detail. We 
regard this solution as a background spacetime for the perturbation analysis. The situation 
can be understood from Fig.l. The "density" D is assumed to have compact support, 
< f < f„, which is depicted by the shaded region. We flnd from Eqs. ()3-3|1 and ()3-7|l that if 
D does not vanish at the symmetric axis r = 0, the components of the stress-energy tensor 
T^y with respect to the coordinate basis diverge there, and the same is true for the Ricci 
tensor by Einstein equations. This is a naked singularity which is depicted by the dashed 
line at r = on the interval < t < in Fig.l. Although all the scalar polynomials of 
the Riemann tensor vanish there, freely falling observers experience an infinite tidal force at 
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r = for < t < fw^^'' The short dashed hne corresponding to t = r represents the Cauchy 
horizon associated with this naked singularity. The region satisfying t > at r = is a 
conical singularity which is depicted by the dotted line in Fig.l. 



Conical 



Singularity 




Intermediate 
Singularity 







Cauchy horizon 



Collapsing Null Dust 



Fig. 1. Morgan's cylindrical null dust solution. Null dust exists in the shaded region. In this 
situation depicted, D(v) > ioi < v < Vy^. The dashed line on the t-axis corresponds to the 
intermediate singularity at which an observer experiences an infinite tidal force, although the 
scalar polynomials of the Riemann tensor do not vanish. The dotted line on t-axis represents the 
conical singularity. The Cauchy horizon is represented by the short dashed line corresponding 
to t = r. 

Now let us consider a general perfect fluid with very large collapsing velocity. To carry 
out the linear perturbation analysis, we introduce a small parameter e and assume the orders 
V — 0{e) and -0 = 0(e). Further, we rewrite the variables 7, it! and D as 



D^DBil + Sn), 
and assume that 5^, 5r and 5d are 0{e), where 

D - ^ 

We call the perturbative analysis with respect to this small parameter e the 'high-speed 
approximation'. 



(3-9) 
(3-10) 

(3-11) 



(342) 
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The first-order equations with respect to e are given as follows: the Einstein equations 
lead to 



5^' = SnGDee'^ {6^-iIj + 6d- 2d^{r6R) - 2PV/Db} + (r^R)", (3-13) 

6^ = SnGDse'^ {6^ - ^ + 5d - 29,(r5^) - V} + (r^^)', (3-14) 

-S-5'' = -^^e^^PV, (3-15) 
r 

rh - [rdR)" = IQnGe^^ {Db - 2P)V, (3-16) 

^-ij" =^e^^DB-2P)V; (3-17) 



the conservation law (j2-16|) leads to 



aASo + s-,-^)--^^-^v--[v -—V 



+ J_|(P^)-_P^;^}; (3.18) 



the Euler equation ()2-17|) becomes 



PV 

du{iDB-2P)V} = —, (3-19) 

where we have used Eq. fl3-18|) . 

The C-energy E up to first order in e is given by 

^ = ^ [1 - e~2^^ + 2e-2T^ {6^ - (r(5«)'}] . (3-20) 

o 

We can easily see in Eq. ()3-12|) that •jb is constant in the vacuum region, Db = 0. Further, 
from Eqs. ()3-13|) and ()3-14|) . we find that 6^ — (r5/j)' is constant in the vacuum region. There- 
fore, up to first order, the C-energy is constant in the vacuum region. This means that in 
the vacuum region, the C-energy flux vector J'^ vanishes up to this order, and thus it is a 
second-order quantity given by 

Tt r f:,.2 , .a 



^j-r = -^i'i''- (3-22) 

In this case, the C-energy flux vector takes a form very similar to that of the massless 
Klein- Gordon field. 



8 



§4. Does pressure prevents the high-speed coUapse? 



Assuming that both the energy density p and pressure p of the fluid are positive, we 
introduce the quantity 

c,{u,v) = ^^ (4-1) 

If the dominant energy condition is satisfled, Cs is smaher than or equal to the speed of 
light. However, in this paper, we assume a more stringent constraint on Cs, namely 

Cs < 1. (4-2) 

As shown below, the present approximation scheme cannot treat the case of Cg = 1. 
From Eqs. ()2-l()j) and ()2-lH) and by using Eq. ()4-H) . we obtain 

P-J^.D. (4.3) 

With the above equation, the Euler equation ()3-19p can be rewritten as 

Vl + Cs^ J (l-Cs^)r 
Integrating this equation formally, we obtain 



T/ 1 + 



dx- 



Juiv) {I - Cs'^ix,v)}{v - x) 



(4-5) 



where U{v) is an arbitrary function of the advanced time v. In the above equation, the 
velocity perturbation V is not deflned for Cg = 1, and thus the high-speed approximation is 
not applicable to this case. 

We consider two cases, one in which Cg is bounded below by some positive constant, and 
one in which Cg vanishes in the large energy density limit, p ^ oo. 

4.1. Bounded case 

In case that Cg is bounded below by some positive value, we have the inequality 

2r 2 

> b\ (4-6) 



1-Cg2 

where b is some non- vanishing constant. Then from Eq. ()4-4p . we obtain 



d^\n[-—^V]> . (4-7) 

1 + Cg / V — u 



9 



Integrating the above inequality, we obtain 

V > Const X r-^\ (4-8) 

where 'Const' should be positive. We find in the above inequality that the velocity perturba- 
tion V diverges when the fiuid element approaches the symmetric axis, r = 0. This implies 
that in this case, the high-speed collapse is necessarily decelerated so significantly that the 
high-speed approximation breaks down. 

4.2. Vanishing case 

If Cs vanishes in the limit of a high energy density, p —>■ oo, the behavior of the perfect fluid 
can differ from the case of bounded Cg. When the fluid elements approach the background 
singularity, i.e., r = in the region of non- vanishing Db, the energy density will become 
larger and larger. Thus, in this case, the asymptotic behavior of Cg near the background 
singularity will be given by 

c,\u,v) ~ a\ty = - ry ~ a\vy, (4-9) 

where Cg is some function, and g is a positive constant. Substituting the above asymptotic 
form of Cg into Eq. ()4-5|) . we obtain 

V ~ Cv^{v) exp l'^Cs\vy\ , (4-10) 



where 

Cv[v) = exp 



{v - U{v)Y 



(4-11) 



Note that Cy should have the same support as the background density variable Db- There- 
fore, in this case, the velocity perturbation remains flnite, even at the background singularity. 
This seems to imply that the naked singularity forms at the background singularity through 
this gravitational collapse. However, the full-order analysis is necessary to reach a deflnite 
conclusion. 

Here, conversely, we consider the asymptotic equation of state that realizes the asymptotic 
behavior ()4-9p . We flnd from Eq. ()2-lU|) that as the background singularity is approached, 
the energy density behaves as 

p ^ = (function of vf x -, (4-12) 

where we have used the fact that V becomes a function of the advanced time f at r = 0, as 
seen from Eq. (j4-l(J|) . Therefore together with Eqs. (j4-H) and (j4-9j) . the above equation leads 
to 

p = ~ Const X p^-^ (4-13) 
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near r = 0, where 'Const' should be positive. This resuh suggests that a very soft equation 
of state is necessary so that the high-speed collapse is not halted by the pressure before 
spacetime singularity formation. 

4.3. Generation of gravitational waves 

Equation ()3-17|) implies that a larger velocity perturbation V leads to a larger \ip\. Be- 
cause the C-energy flux in the asymptotically flat region is given by Eq. (j3-22| ). this suggests 
a large gravitational emission by large V. 

In the bounded Cs case, the growth of the velocity perturbation V is unbounded as the 
background singularity is approached. Although a velocity perturbation V larger than unity 
is meaningless in the present approximation scheme, this result suggests that a significant 
amount of gravitational radiation is generated by pressure deceleration near the background 
singularity. This is consistent with Piran's numerical results, because Piran's numerical 
simulations for a cylindrically distributed ideal gas show that a large amount of gravitational 
radiation is emitted when the pressure bounce occurs. However, here we stress that in this 
paper, we study the result for an arbitrarily large initial collapsing velocity of the fluid, while 
Piran's numerical study is restricted to cases with momentarily static initial data. 

In the case of vanishing Cg, the velocity perturbation takes a finite value even at the 
background singularity. Therefore the situation here is similar to that for a dust fluid in the 
generation processes of gravitational waves. 



When the gravitational collapse realizes a very high energy density, the fluid will behave 
as a mono-atomic ideal gas, due to the asymptotic freedom of the elementary interactions. 
Furthermore, the equation of state will be similar to that of the radiation, i.e., with p = 
p/3. This corresponds to the bounded Cg case, and the high-speed collapse is necessarily 
decelerated so significantly that the high-speed approximation becomes inapplicable before 
the fluid elements reach the background singularity. This result seems to imply that the 
formation of spacetime singularities is prevented by the occurrence of the pressure bounce. 
However, even if this is the mentioned in §1, it is physically important in a practical 

sense only whether the spacetime border forms before the occurrence of the pressure bounce. 
Therefore, here we consider the formation of the spacetime border in the bounded Cg case. 

We consider a linear equation of state, i.e., Cg =constant, which includes the case of 
radiation, Cg^ = 1/3. Using Eq. ()4-3p . the Euler equation ()3-19|) becomes 



§5. Appearance of the spacetime border 



duV 



2cgV 



(5-1) 
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As mentioned previously, the above equation is singular for stiff matter with Cg = 1. The 
case Cg = 1 corresponds to a massless Klein-Gordon field, and it is easy to see that no 
gravitational radiation is generated in the cylindrically symmetric case. 
Equation (j5-l|) can be easily integrated, and we obtain 

V = CL{v)r~^''''^^-''\ (5-2) 

where Cl{v) is an arbitrary positive function of the advanced time f, and it should have 
the same support as the background density variable Db- As discussed in the preceding 
section, V diverges to +oo at the background singularity in the case that < Cg < 1. This 
implies that the perturbation analysis necessarily breaks down in the neighborhood of the 
background singularity. 

Once we know the velocity perturbation V , we can easily obtain solutions for 5^, 
and Tp from Eqs. (j3-13|) — (I3-I7p using the ordinary procedure to construct solutions of wave 
equations in flat spacetimes. We can also easily obtain a solution for do by solving Eq. ()3-18p . 
However, in this paper, we do not construct the solutions for these variables but instead focus 
on the velocity perturbation V only. 

The high-speed approximation scheme is applicable only to situations in which V <^ 
1. Therefore, Eq. ()5-2|) leads to a necessary condition for the applicability of high-speed 
approximation, 

r > Cl^^-'"^'^'^\ (5-3) 

where r is the radial position of a fluid element. It should be noted that if the arbitrary 
function Cl is chosen sufficiently small, the present approximation scheme may be applicable 
even in the case of very small r. Because the tidal force experienced by freely falling observers 
diverges as the background singularity at r = is approached, ^^-^ the velocity perturbation 
V remains sufficiently small for the case of an arbitrarily large tidal force for freely falling 
observers ii Cl is chosen sufficiently small. This implies that if the initial collapsing velocity 
is sufficiently large, the gravitational collapse of a cylindrically distributed perfect fluid with 
a linear equation of state forms a spacetime border. Because there is no event horizon in the 
present case, the spacetime border is necessarily visible. This result implies practically the 
violation of cosmic censorship. 

§6. Summary and discussion 

In this paper, we have presented a high-speed approximation scheme applicable to a 
cylindrically symmetric perfect fluid. Within this approximation scheme, we have found 
that in the case of Cg bounded below by some positive value, the high-speed collapse of the 
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perfect fluid is necessarily decelerated so significantly that the high-speed approximation 
becomes inapplicable. This result reveals the possibility of pressure bounces, because the 
deceleration of the collapsing velocity implies that the pressure repulsion overwhelms the 
gravitational attraction. If this expectation holds, then the spacetime singularity formation 
resulting from the adiabatic collapse of a cylindrical perfect fluid of a physically reasonable 
ideal gas is impossible. However, here we should note that if the initial collapsing velocity 
is very large, a spacetime border (i.e., a region with spacetime curvature so large that the 
quantum effects in gravity are important) is realized before the high-speed approximation 
scheme breaks down. Because a horizon does not form in the system considered here, the 
spacetime border is necessarily visible. Therefore, this result implies practically the violation 
of cosmic censorship. By contrast, in the case that Cg vanishes in the limit of high energy 
density, p ^ oo, the high-speed collapse continues until a naked singularity forms. Thus, 
it is likely that if the equation of state is sufficiently soft, so that p is proportional to p^^^ 
with < g < 1, high-speed collapse is not prevented by the effect of the pressure, and thus 
a naked singularity forms. 

Equation ()3-17|) implies that a large velocity perturbation V leads to the efficient genera- 
tion of gravitational radiation. Therefore, the unlimited growth of V in the bounded Cg case 
implies the generation of significant gravitational radiation. Contrastingly, in the case that 
Cg vanishes in the high energy density limit, V is finite, even at the background singularity. 
In this case, the behavior of the perfect fiuid will be similar to that of a dust fiuid, in which 
case a thinner cylindrical shell generates the larger amount of gravitational radiation. 

In the case that Cg vanishes in the high energy density limit, the gravitational emission 
from the naked singularity might have some physical importance. In the present approxi- 
mation scheme, the metric variables are regarded as test fields in the background spacetime. 
Thus, if we appropriately impose boundary conditions on these variables at the background 
singularity, we can obtain solutions of the linearized Einstein equations and estimate the 
energy carried away by the gravitational waves from the naked singularity in this case. This 
topics will be discussed elsewhere. '^^^ 

The approximation scheme used here is a kind of perturbation method. Strictly speaking, 
the full-order analysis is necessary to reach a definite conclusion, although the results up to 
first order might have significant physical meaning. Therefore, a higher-order analysis should 
be meaningful, but this is a future work. 



13 



Acknowledgements 



We are grateful to our colleagues in the astrophysics and gravity group of Osaka City 
University for helpful discussions and criticism. KN thanks A. Hosoya and A. Ishibashi 
for useful discussion. This work is supported by a Grant-in-Aid for Scientific Research 
(No. 16540264) from JSPS and also by a Grant-in-Aid for the 21st Century COE "Center 
for Diversity and Universality in Physics" , from the Ministry of Education, Culture, Sports, 
Science and Technology (MEXT) of Japan. 



References 



C. Cutler and K. S. Thorne, in Proceedings of the 16th International Conference on 
General Relativity and Gravitation, edited by N.T. Bishop and S.D. Maharaj (World 
Scientific, 2002), p. 72. 

H. Tagoshi et al., Phys. Rev D63 (2001), 062001. 
M. Ando et al., Phys. Rev. Lett. 87 (2001), 3950. 

T. Nakamura, M. Shibata and K. Nakao, Prog, of Theor. Phys. 89 (1993), 821. 
R. Penrose, Riv. Nuovo Cim. 1 (1969), 252. 

H. Iguchi, T. Harada and K. Nakao, Phys. Rev. D57 (1998), 7262. 

H. Iguchi, T. Harada and K. Nakao, Prog. Theor. Phys. 101 (1999), 1235. 

H. Iguchi, T. Harada and K. Nakao, Prog. Theor. Phys. 103 (2000), 53. 
K. Nakao, H. Iguchi and T. Harada, Phys. Rev. D63 (2001), 084003. 

T. Harada, H. Iguchi and K. Nakao, Prog. Theor. Phys. 107 (2002), 449. 

T. Harada and K. Nakao, Phys. Rev. D70 (2004), 041501. 

N. Arkani-Hamcd, S. Dimopoulos and G. Dvali, Phys. Lett. B429 (1998), 263. 

I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B436 
(1998), 257. 

L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999), 3370. 

L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999), 4690. 

P. Joshi, Global Aspects in Gravitation and Cosmology, (Oxford University Press, 

1993). 

K.S. Thorne, in Magic Without Magic; John Archibald Wheeler, edited by J. Klauder 

(Prieman, San Prancisco, 1972), p231. 

S. A. Hayward, Class. Quantum Grav. 17 (2000), 1749. 

T. Piran, Phys. Rev. Lett. 41 (1978), 1085. 

F. Echeverria, Phys. Rev. D47 (1993), 2271. 



14 



T. Chiba, Prog. Theor. Phys. 95 (1996), 321. 

K. Nakao and Y. Morisawa, Class. Quantum Grav. 21 (2004), 2101. 

T. A. Morgan, Gen. Rclat. Gravit. 4 (1973) 273. 

T. Apostolatos and K.S. Thorne, Phys. Rev. D46 (1992), 2435. 

P. S. Letelier and A. Wang, Phys. Rev. D49 (1994), 5105. 

B.C. Nolan, Phys. Rev. D65 (2002), 104006. 

P. R. C. T. Pereka and A. Wang, Phys. Rev. D62 (2000), 124001. 
B.C. Nolan and L.V. Nolan, Calss. Quantum Grav. 21 (2004), 3693. 
K.S. Thorne, Phys. Rev. 138 (1965), B251. 

M.A. Melvin, Phys. Lett. 8 (1964), 65; Phys. Rev. 139 (1965), B225. 
R. M. Wald, General Relativity (The University of Chicago Press, 1984). 
K. Nakao and Y. Morisawa, in preparation. 



15 



